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1. Introduction 

A topological charge is one of the most fundamental quantity in QCD. It characterizes the 
vacuum structure. Lattice QCD is a main tool in the study of the topological charge Lattice 
QCD allows us to perform a nonperturbative analysis in a systematic way. 

The topological charge is often measured with a gluonic field strength operator on the lattice. 
Though it suffers from noisy ultraviolet fluctuations, a smoothing technique tames them so that a 
discernible signal of the topological charge can be obtained. Cooling or smearing have been used 
for smoothing gauge fields. Recently, a gradient flow is also employed. In contrast to the traditional 
cooling and smearing, the gradient flow has an advantage that it provides a continuous change of 
the gauge field. The gradient flow accomplishes a better control of smoothing. 

Alternatively, the topological charge can be calculated with a fermionic definition. The topo¬ 
logical charge is determined, for example, by the index theorem with the overlap-Dirac operator. A 
clear advantage of the fermionic definition is that the result is guaranteed to be an integer. A subtle 
point, on the other hand, is the integer values depend on the choice of the definition, due to a finite 
lattice spacing. In the case of the overlap-Dirac operator, a value of the topological charge occa¬ 
sionally changes according to a parameter in the formulation. Consistency check of the topological 
charges in the fermionic and gluonic definitions would be helpful as an estimator of the scaling 
violation. 

In this work, topological charges are computed with gluonic operators on Nf = 2 topology 
fixed gauge configurations. The measurements are performed using several smoothing techniques. 
Cooling with plaquette and improved local actions, APE and HYP smearing, as well as gradient 
flows are employed. The results are compared with each other, and with the values obtained using 
the overlap-Dirac operator. Similar attempts are reported in Refs. |^. 

2. Setup 

2.1 Gauge configuration 

Measurement of the topological charge is performed on Nf = 2 gauge configurations provided 
by JLQCD Collaboration [^. The lattice size is 16^ x 32 at the lattice spacing of a = 0.118(2) fm. 
The gluon action is Iwasaki-type improved gauge action. 



(2.1) 


where j3 = 6/gQ, Cq = 3.648, cf = —0.331. The quark action is an overlap-Dirac fermion action. 



where m is the bare quark mass. Dw(—mo) is the Wilson operator with a negative mass, —rtiQ = 
— 1.6. Furthermore, unphysical Wilson fermion y/ with a negative mass as well as twisted mass 
terms are added to fix the topological charge defined by the index theorem 2index> 


5Sw = VDw{—mo)V + + inY5‘^3)^ 


(2.4) 
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j3 m 

Gindex 

# conf 

MD time 

2.30 0.05 

-2 

50 

2500 


Table 1: Parameters of the gauge configurations. Molecular Dynamics time is the number of trajectories 
multiplied by the trajectory length. 


where 0 is a pseudofermion, and /r is the twisted mass parameter, /r = 0.2 is employed in the 
configuration generation. Each configuration is separated by 100 trajectories with its trajectory 
length 0.5. The main simulation parameters are summarized in Table ||. 


2.2 Gluonic topological charge operator 


The topological charge is measured using a gluonic field strength. 


Gimprove = , 

(2.5) 

/ilxl,2 _ ^ V'’ ^ p.1 X 1,2 p-1 X 1,2 

^ “ 219^2 ^pG 5 

p,V,p,G 

(2.6) 

pi X 1,2 _ ^ rptlxl,2iAH prAH _ ^ /pi ^ 

(2.7) 

where is the cloverleaf constructed with 1x1 plaquettes, and with 1x2 rectangular 

loops. The improvement coefficients cq and ci can be tuned to reduce the scaling violation in the 
topological charge operator. Three types of (co,ci) are investigated. 

(co,ci) = (1,0) Naive-type, 

(2.8) 

= (5/3,—1/12) Symanzik-type, 

(2.9) 

= (3.648,-0.331) Iwasaki-type. 

(2.10) 


Figure |i| displays ci dependence of the topological charge on a single configuration smoothed by 
Wilson flow. The topological charge wifh Symanzik-fype coefficienfs has fhe smallesf deviafion 
from an integer. Since fhe deviafion is originated from fhe finife laffice spacing, if implies an 
efficienl reducfion of fhe scaling violafion by Symanzik-fype operator. Calculafions on ofher con- 
figurafions show a similar fendency. Based on fhis resulf, Symanzik-fype coefficienfs are employed 
in fhis work. 


2.3 Smoothing 


Three kinds of smoothing techniques are evaluated: cooling, smearing, and gradient flow. 
Smoothing is required to suppress noisy ultraviolet fluctuations, while keeping a topological struc¬ 
ture. Although any smoothing is expected to give a consistent result in the continuum limit, it is 
valuable to find a method that has the least lattice artifact. 

Cooling eliminates ultraviolet noises by replacing each link variable such that the local action 
is minimized [Q]. For the local action, not only a naive plaquette action, but also Symanzik and 
Iwasaki actions are employed with the coefficients of Eq. (|2.8[)-(Tl^. 
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Figure 1: Improvement coefficient ci dependence of improved topological charge ^improve on a single gauge 
configuration. Difference of ^improve from an integer is plotted. 


Another way of smoothing is smearing. APE smearing is defined by 






P™jsu(3) ~ +—L^{x) 

^ Uyix)U^{x + v)Ul{x + ^) 


,a = 0 . 6 . 


±v/p 

In addition, HYP smearing is also examined. 


( 2 . 11 ) 

( 2 . 12 ) 


ujclix) 

a\ 


P™jsu(3) 

Projsu(3) 

P™jsU(3) 


(l-ai)f/^(x) + ^ £ U^l{x)u''^.l{x + v)ul^}f{x + fi) 




^(2).t/ 
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(l-a2)f/pW + ^ L ufiv{x)U^Xi^ + P)Up}lli^ + ^) 

^ ±p^p,v 


(1 — a3)17^(x) + ^ ^ Ua{x)U^{x +a)ul{x +n) 


±afp,v,p 


0.75, a2 =0.6,a3 =0.3. 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


Smeared gauge fields are projected back to SU (3) by Maximum SU (3) projection. 

Proj“f;'='(C/„{x)) = Re Tr (2.17) 

An alternative smoothing is given by the gradient flow |^. The evolution of the gauge field is 
determined by 


V^{x,t=0) = U^{x), 


(2.18) 

(2.19) 


where t is the flow time, and S is an action without its coupling constant. Similar to the cooling 
case, plaquette, Symanzik, and Iwasaki actions are employed. The flow equation is solved by the 
fourth order Runge-Kutta in the commutator-free method 0. The Runge-Kutta step size dt is 
chosen to be 0.02. The systematic error associated with discretization of the flow time is definitely 
below the statistical error. 
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Figure 2: Smoothing step dependence of ^improve using cooling (top panels), smearing (middle panels), and 
gradient flow (bottom panels). 


2.4 Results 

Figure ^ illustrates cooling and smearing step dependence of the improved topological charge. 
The flow time dependence is also plotted. The flow time is multiplied by a factor of three, which 
is expected from a perturbative analysis [^. In every case, the topological charge has an integer 
value with a sufficiently large number of steps. A small number of smoothing steps leads to a 
fake plateau i.e. a semi-stable value of the topological charge. The number of smoothing steps is 
determined to satisfy the admissibility condition, max[ReTr(l — t/piaq)] < 0.067 |[^, |^]. In Fig. ||, 
Wilson flow time dependence of the plaquette is shown. No jump of the topological charge seem to 
be triggered, if the admissibility condition is fulfilled. It should be mentioned max[ReTr(l — t/piaq)] 
does not always decrease as the flow time grows, though the value summed over the spacetime falls 
off monotonically. 

Figure ^presents histograms of the improved topological charges. Since the topological chai'ge 
determined by the index Qindex is fixed in the configuration generations, the histogram is expected to 
have a sharp peak around Qmdes, supposing the scaling violation is small. The histograms obtained 
by cooling with improved local actions show the expected behavior. Almost all of the topological 
charges agree with Qindex- On the other hand, cooling using the plaquette action has a broad his¬ 
togram. It implies a relatively large lattice artifact in the unimproved cooling method. Analogous 
trends are observed in other smoothing procedures. HYP smearing has a narrow histogram, while 
APE smearing does not. Symanzik and Iwasaki flows form a sharp peak in the histogram. On the 
contrary, Wilson flow brings a wide peak. Improved smoothing methods leads to higher consis- 
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Figure 3: Wilson flow time dependence of maximum of (1 — plaquette) (left panel), and the value summed 
over the spacetime volume (right panel). 


tency with Qindex; indicating the scaling violation in the topological charge is suppressed well by 
the improvement. 

2.5 Conclusion 

Systematic comparison of topological charges is presented. Topological charges are measured 
on Nf = 2 topology fixed configurations. Several smoothing techniques are evaluated using a 
gluonic topological charge operator of Symanzik-type coefficients, which give a topological charge 
with the smallest deviation from an integer. 

Cooling with improved actions, HYP smearing, and improved gradient flows are found fo be 
advanfageous. More fhan 90% of fhe topological charges are consisfenf wifh fhose obfained by fhe 
index fheorem. If indicates fheir lattice arfifacfs are reduced efficienfly. On fhe ofher hand, cooling 
wifh plaquette acfion, APE smearing, and Wilson flow lead fo parfial mafches. The agreemenf is 
limifed fo 70-80%. Scaling violations seem fo be comparatively large in fhese smoofhing mefhods. 

Scaling properties as well as finife size effecfs of fhe fopological charge have nof been inves- 
figafed. If is imporfanf fo esfimafe fhem, buf is beyond fhe scope of fhis work due fo limifafion of 
fhe gauge configurations. The gauge configurafions have been generated af a single lattice spacing 
and spafial volume. These evaluations are lefl for fhe fulure work. 


Acknowledgments 


I would like fo fhank members of MEXT SPIRE Eield 5 Projecf 1 and Bridge-i-i- develop- 
menf learn, as well as H.Eukaya and K-I.Nagai for valuable discussions. I am grateful to JEQCD 
Collaborafion for providing fheir gauge configurations [^. This work is in pari based on Bridge-i-i- 
code [12, lA, ^] . This work is supported by JEDG conslrucled over SINET3 of Nil, MEXT SPIRE 
and JICEuS, and Granls-in-Aid for Scientific Research Granl Number 24540250. 


References 

[1] For a recent review, see M. Mueller-Preussker, PoS (LATTICE2014), 003 (2014). 


6 























Comparative study of topological charge 


Yusuke Namekawa 


50 

45 

40 

35 


20 

15 

10 


50 

45 

40 

35 


20 


16^x32 

Coolinq(plaq) ■ 


16^x32 

Coolinq(SYmanzik) B 

45 

16^x32 

Cooiinq(lwasaki) B 

^=2.30, mud = 0.050 


2.30, mud = 0.050 


^=2.30, mud = 0.050 


Qindex =-2 

40 

Qindex =-2 


40 

Qindex =-2 




35 



35 





30 



30 





20 



20 





15 



15 




1 

10 



10 




1 - 



■ 

5 



5 -5 -4 -3 - 2-10 1 


-5 -4 -3 

2-101 


-5 -4 -3 - 

2-101 




16“x32 


16^x32 

Symanzik flow b 


16^x32 

Iwasaki flow b 

2.30, mud = 0.050 


2.30, mud = 0.050 



2.30, mud = 0.050 


Qindex =-2 

40 

Qindex =-2 


40 

Qindex =-2 



35 



35 




30 



30 




20 



20 




10 


_■_^^_ 

10 

5 
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